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Abstract: We present a next-to- leading order calculation of heavy flavour production in 
hadronic collisions that can be interfaced to shower Monte Carlo programs. The calculation 
is performed in the context of the POWHEG method [1]. It is suitable for the computation 
of charm, bottom and top hadroproduction. In the case of top production, spin correlations 
in the decay products are taken into account. 
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1. Introduction 

The study of heavy flavoured hadronic final states produced in hadron collisions has re- 
ceived great attention in the past twenty years. This interest was driven partly by the 
search for the top quark, and partly by the circumstance that processes characterized by 
one single large energy scale (the heavy quark mass in this case) are reasonably controlled 
by perturbative QCD due to asymptotic freedom. Theoretical progress has therefore par- 
allelled the experimental efforts in this direction, and has by now reached a remarkable 
level of refinement. 

The next-to-leading order (NLO) calculation of heavy flavour hadroproduction has 
been matched to the HERWIG SMC generator [2, 3] in the context of the MC@NLO 
formalism [4, 5]. This development is particularly important, since it has always proved 
difficult to correctly simulate heavy flavour production in standard shower Monte Carlo 
programs (SMC from now on), because of the impact of higher order processes like flavour 
excitation and gluon splitting. 

In ref. [1] a method for interfacing NLO calculations with SMC generators was sug- 
gested, that overcomes some drawbacks of the MC@NLO technique. In particular, the 
implementation of a given production process with the method of ref. [1] is independent of 
the SMC it will be matched to, whereas this is not the case of MC@NLO, which requires an 
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SMC-dependent (but process independent) contribution which is very laborious to obtain 
(this is the reason why so far MC@NLO has been matched only to HERWIG). On the 
other hand, the technique of ref. [1] requires the use of a soft shower (the vetoed-truncated 
shower) which is not available in present SMC generators, but whose effects are expected 
to be small. Furthermore, the approach of ref. [1] generates events with positive weights, 
while MC@NLO generates also events with negative weight. For this reason the method 
has been referred to as POWHEG, for POsitive Weight Hard Event Generator. In ref. [6] 
the POWHEG method was successfully employed in the calculation of Z pair production 
in hadronic collisions. 

In the present work, we apply the POWHEG method to heavy flavour hadroproduction. 
NLO cross section formulae for this process have been available for a long time [7, 8, 9]. In 
the application of the POWHEG method, no new analytic results are needed. Furthermore, 
since the kinematics of the process is similar to that of Z pair production [6] , many results 
of that paper directly extend to the case at hand. The result of our work provides a 
valuable alternative to the MC@NLO program. It can be used with both HERWIG and 
PYTHIA [10], and it should not be difficult to interface it to other SMC's as well. It 
furthermore provides a viable alternative, in cases when negatively weighted events in 
MC@NLO cause severe problems, like, for example, when the heavy quark is relatively 
light compared to the energies involved. 

The web location of the code repository can be found in the manual, ref. [11]. 

2. Kinematics and cross section 

The differential cross section for the production of heavy quark pairs in hadronic collisions 
was computed in ref. [7], [8] and [9] up to order o? s . In this section, we formulate the result 
of ref. [9] in a form which is suitable for the generation of the hardest emission using the 
procedure proposed in ref. [1]. 

The order-a^ cross section for the process H1H2 — > QQ + X can be written as the sum 
of four terms: 

da = daV> + da {s ^ + da® + da {c \ (2.1) 

Here da^ is the leading-order (Born) cross section. The term da^ collects order-aj? 
contributions with the same two-body kinematics as the Born term, namely one-loop cor- 
rections and real-emission contributions in the soft limit. Finally, da^ represents the cross 
section for real emission, and the corresponding subtractions, in a generic configuration, 
while da^ is a remnant of the subtraction of initial-state collinear singularities. 

2.1 Born and Soft- Virtual kinematics 

At leading order, the relevant parton subprocesses are 

g(pi)+g(P2) Q(ki) + Q(k 2 ), 

q(pi) + q(P2) - Q(h) + Q(k 2 ), (2.2) 
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where q is a quark or antiquark of any flavour, q the corresponding antiparticle, and 
g is a gluon. Particle four-momenta are displayed in brackets; we have p\ = p\ = 0, 
k\ = k 2 = m 2 Q , where m Q is the heavy quark mass. 

The kinematics of heavy quark pair production is entirely analogous to that of vector 
boson pair production. We introduce the reduced Mandelstam invariants 1 

s = ( Pl + p 2 ) 2 , t = (pi - h) 2 - ml, u = (pi - k 2 ) 2 - m 2 Q , (2.3) 

related by s + t + u = 0. 

Event generation is conveniently performed in terms of the invariant mass M q q and 
the rapidity Y q q of the QQ pair in the laboratory frame. They are given by 

M 2 Q = (k 1 + k 2 ) 2 = Xl x 2 S (2.4) 

v 1, (gi + P2) + (gi + P2) 3 1, x 1 
2 (pi +p 2 r ~ (Pi +P2) 6 2 x 2 

where S is the squared center-of-mass energy of the colliding hadrons, and xi, x 2 are the 
fractions of longitudinal momenta carried by the incoming partons. Equations (2.4,2.5) 
yield 



M 2 - M 2 - 1 

xi = W e Y QQ = x b i; x 2 = W e- y «<5 = x 62 ; dxidx 2 = - dY QQ dM 2 Q . (2.6) 

We adopt as two-body kinematic variables the set v = {Mqq, Y Q q, cos #1} (which we will 
call the Born variables henceforth), where 9i is the angle between pi and ki in the partonic 
center-of-mass frame, so that 

M 2 n 

t = ^(1- /?cos0i), (2.7) 



where 



4m? 



/3=7W; P=Tj^- (2-8) 

QQ 

Using eq. (2.6) and the usual expression of the two-body phase space measure d& 2 , one 
immediately finds 

d§ 2 dxi dx 2 = d cos 61 dM 2 Q dY QQ . (2.9) 



In order to keep our notation similar to that of ref. [1] , we define 

d<$> v = dcosOi dM 2 Q dY QQ . (2.10) 
The appropriate integration region for the variables v is 

4m 2 Q<M 2 Q<S, -log-|« <Y Q(3 <--log-p, -1<cos0i<1. (2.11) 



1 We depart slightly from the notation of ref. [6], where t and u have different definitions. 
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The Born cross section is given by 



da {h) = d$ v Y^Bj{v,n), (2.12) 

3 

where 

Bj(v,n) = ^f^{x hl ^)f^{x b2 ,^)Mf_ j {M 2 QQ ,t). (2.13) 

The index j represents the light quarks and antiquarks (j / 0) and the gluon (j = 0). It 
ranges between — ne and ri£, where ri£ is the number of light flavours. In the following we 
will also use q to represent all values of j except for j = 0, q to represent — j, and g to 
represent j = 0. f!?(x,ii) denotes the distribution function of parton j in the hadron H, 
and fi is a factorization scale. The function A4^}_j(s,t) is the squared invariant amplitude, 
summed over final-state polarizations and averaged over initial-state polarizations and 
colours, divided by the relevant flux factor. We have 

where N c = 3, T P = 1/2, C P = 4/3, £> A = 8. 

Order-a^ contributions to the cross section arise from one-loop corrections to the two- 
body process eq. (2.2), and from real-emission subprocesses at tree level. The contribution 
of one-loop diagrams must be summed to the one-gluon emission cross section in the soft 
limit, in order to obtain an infrared-finite result. The resulting contribution, usually called 
the soft-virtual contribution, has the same kinematic structure as the leading-order term: 

da^ = d$ v ^Vj(v,n), (2-15) 

3 

where 

Vj(v,») = j^f? 1 (x bl ,rif^(x b2 ,t,)M^.(MZ Q ,t,t, 2 ). (2.16) 

The invariant amplitude A^ s ^-(M|g, t, fi 2 ) is the sum of the virtual corrections to the 2 
body subprocesses, and the soft contribution of real emission. The sum is finite, and it was 
computed in ref. [9]. 

2.2 Real emission kinematics 

We consider now the real-emission subprocesses 

q(pi) + q(P2) - Q(fci) + Q(k 2 ) + g(k) (2.17) 

q(pi) + g(P2) - Q(fci) + Q(k 2 ) + q(k) (2.18) 
g(pi) + ?(P2) ^ Q(fci) + Q(k 2 ) + (2.19) 
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in a generic kinematical configuration. The processes (2.17-2.19) are characterized by five 
independent scalar quantities, which we choose to be 



s = (pi+ P2) 2 , t k = (pi - k) 2 , u k = (p 2 - k) 2 , 
qi = (pi - h) 2 - m 2 Q , q 2 = (p 2 - k 2 ) 2 - m 2 Q , 



as in ref. [9]. We introduce the variables 



x 



y = cos 9, 



(2.20) 
(2.21) 

(2.22) 



where 9 is the scattering angle of the emitted parton in the partonic center-of-mass system. 
With these definitions, 



t k = --(!- x){l-y)- 



u k = - gt 1 - x )( 1 + y)- 



(2.23) 



It is easy to show that in the case of the subprocesses (2.17-2.19) one has 



and therefore 



Xl 



Xl s + u k 


= ^log 


x 2 s + t k _ 





ai2- (l-x)(l + y) 
x 2 2-(l-x)(l-y) 



M, 



QQ 



x x\x 2 S, (2.24) 



Xbi 2-(l-x)(l-y) _ 

y/E\ 2-(l-x)(l+yY 



x 2 



Xb2 



/2-(l-x)(l + y) 



v ^V2-(l-x)(l-y) 



and 



dx x dx 2 = — dM 2 QQ dY QQ . 



(2.25) 



(2.26) 



The range for the variable x is restricted by the requirement that both x\ and x 2 be smaller 
than one; this gives 

Xmin < x < 1, (2.27) 

with 



max 



2(l+!/)a& 



^1 + x 2 bl ) 2 (l - y) 2 + 16yx 2 bl + (1 - y)(l - x 2 hl ) 

2(l-j/)xg 2 

yj(l + x 2 2 ) 2 (l + yf - lSyx 2 2 + (1 + y)(l - x 2 



(2.28) 



Note that x m i n depends explicitly on y, and implicitly on Mqq and Y Q q through Xbi,Xb 2 . 
It can be checked that x m [ n is always larger that Mgg/S, as required by the definition of x. 

In the center-of-mass frame of the QQ system, the four-momenta of the produced 
heavy quarks can be parametrized in terms of two angles 9\,9 2 : 

Mqq 



k 2 = ^ 



(1, P sin 9 2 sin 0i, (5 cos 9 2 sin 9\ , [3 cos 9\ 



'1, —(3 sin 9 2 sin 6*i, — (3 cos 9 2 sin#i, — /?cos ^1), 



(2.29) 
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with (3 given in eq. (2.8). Both 9\ and O2 range between and it. Thus, in addition to the 
Born variables v, we have now the three radiation variables r = {x,y, 62}, with 



Zmin <X<1, ~l<y<l, < # 2 < 7T . 

Following ref. [1] we define the corresponding integration measure 

dQ r = dx dy d02 ■ 
From the computation of ref. [9] we obtain 



(2.30) 
(2.31) 



da^ = d<Z> v d$ r I R gg (v, r,fi) + ^2 [ R qq( v > r > A 4 ) + Rqg(v, r, fl) + Rgq(v, r, fl)] \ , (2.32) 



where 
Rgg(v,r, n) 

R q q(v,r,[i) = 

R qg( V i r i 
Rgq(v, r, li) 



(4tt) 2 64tt 2 M% q S \1-x 



+ 



fg 1 ( x l > A*) fg 2 ( x 2 ,li)fgg(x,y,Oi,02,li) 



(4vr) 2 Utt 2 MI q S \\-x 



i-y 



i + y 



(2.33) 



fq 1 («1 , M) /f 2 (»2, /<re(z, V, #1 , #2, A 1 ) 



(4tt) 2 64tt 2 M2 (5 S' Vl-a:/ VlTy/ + 



/«f 1 (zi , m) /«f 2 (#2 , m) / 9S (x,y,0i,9 2 ,fj.) 



cH 2 



(4vr) 2 64vr 2 M 2 (3 5 ^1 - xj \1 T 



fq 1 (.Xl, fi) /f 2 (x 2 , M) fgq(x,y,0i,9 2 ,fi). 



+ 



The functions denote the regularized real emission cross sections for the 

different subprocesses. The functions f qq - and f qg are regular in the limits of soft (x = 1) 
or collinear (y = ±1) emission; they are defined as in eq. (3.3) of ref. [9]. The distributions 
1/(1 — x) p and 1/(1 ± y) + are defined by 



dx g(x) 



1 



1 — x 



1 — x 

%)-Mti) 



(2.34) 
(2.35) 



2.3 Collinear remnants 

The remnants of the collinear subtraction must also be added to get the full cross section. 
This contribution has the form [9] 

dcr {c) = d<5> v dxdy (2.36) 

{ [Lg 9 (v,x,iJ,) {L^q(v,x,n) +L+g{v,x,n) +L+g(v,x,n)) 8(1 - y) 

Q 

+ [ L g~g(v,X,n) ( L qq(v,X,fi) + L~ q -(v,X,n) + L qg (v , X , fl)) 8(1 + y) }, 
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where 



a s (3 
2vr WttS 



J2\ [(l-x)P ii ,(x,0)] 



1 — X 



log- 



M, 



QQ 



+ 2 



log(l — x) 
1 — X 



[(l-x)PlAx,0)] l^gCM^,*)/^^^)/^ (x 2>M ) (2-37) 



and 



log 



/log(l-x)- 



1 — X / Xjl^ 



+ 2 



1 -x 



-[(l-s)J*,(x,0)] (^) [^(M^,*)/^^^)// 2 (x 2 ,/x) (2.38) 



where Pij(x, e) are the leading-order Altarelli-Parisi splitting functions in d = 4 — 2e di- 
mensions, P(j(x,e) their first derivatives with respect to e, and t is given in eq. (2.7). 
From eq. (2.28) we see that the integration range becomes Xbi < x < 1 when y = 1, and 
Xfe2 < a; < 1 for y = — 1. 



3. Cross section for the hardest emission 



The POWHEG method, when applied to a generic process, may require a separated treat- 
ment of each singular region. In the present case (as in the case of Z pair production) 
this is not needed. Our choice of variables v, r is adequate for both collinear regions at the 
same time, the only difference being the sign of y. We have instead to pay attention to the 
flavour structure of the process. In ordinary SMC codes, the flavour structure of the Born 
subprocess is not altered by subsequent radiation. On the other hand, if the hardest radi- 
ation is produced in the context of a NLO calculation, the association of the NLO process 
with a Born subprocess is not always obvious. A given real-emission contribution must be 
associated with its underlying Born process, i.e. the Born process in which it factorizes 
in the collinear limit. In the present case, the collinear regions for the real qq subprocess 
always factorize in terms of the qq underlying Born, and the collinear regions for the real 
gg subprocess factorize in terms of the gg Born process. For the qg (gq) processes, there 
are instead two possibilities (see fig. 1): the underlying Born process is gg (qq) for the 
y = 1 collinear region, and qq (gg) for y = —1. This is the reason why in eq. (2.33) we 
have separated the two collinear contributions and R^. Thus, for a given flavour q, 
we lump together the qq, the qg and the gq real-emission subprocesses. 
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Figure 1: Different underlying Born structure in the two collinear regions of the qg subprocess. 



Following ref. [1], we write the cross section for the event with the hardest emission as 
A g (0) + A q (k T ) ^ r ' ^ + r ' ^ + ^ r ' M d* r 



da = y~]B q (v,iJ, v )d$ v 
i 

+B g (v,/i v )d<f> 



B q (v,n r ) 



Rgg{v, V, Hr) + J2q ( Rtgi v , r , Mr) + Rgq(v ', T, fi, 

Ag(0) + Ag(k T ) . ! '_ d $ r 



B q (v,/J, r ) 



(3.1) 



where Rij is obtained from i?^- by removing the + prescriptions that regularize the x and y 
singularities. The Rij are thus the unregularized real emission cross sections (corresponding 
to R in the notation of ref. [1]). Notice that the real emission contributions having the 
same underlying Born configuration are grouped together in eq. (3.1). Furthermore, 

B q (v,fi) = B q (v,n) + V q (v,fjt) + J d$ r [R qq -(v,r,n) + R qg (v,r,fi) + R+ q (v,r,n)] 
+ dy dx [L+ q (v,x,fj,) + L+g(v,x,fj,)] 5(1 - y) 
+ dy dx [L~g(v,x,fi) + L~ g (v,x,fi)] S(l + y) (3.2) 

J —1 J X m [ n 



B g {v,n) =B g (v,n) + Vg(v,n) + j I d$ 



fl f-1 
+ I dy dx 

-1 J x m - ln 



+ 



f dy f 

J 1 J x m - ln 



dx 



R gg (v, r,n) + Y^ {Rtg( v > r > A*) + R gq( v , r , M)) 
1 

L+ (v,x,fi) +^L+ g (v,x,fi) 6(1 -y) 
g 

L qq(v,X,^) +^2Lg-q(v,X,fl) 5(1 + y) 



(3.3) 



A q (p T ) = exp 



Rqq(v, r,fl r )+ R qg (v, T, /v) + Rgq(v, T, fl 



Bq(v,Ll r ) 



*(v,r) -p T )d$ r 



(3.4) 
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A s(Pt) = exp - 



/ 



Rgg(v, r, Av) + Eg gjgK r , Mr) + E g Rgq(.v, r, Hr) 
Bg(v,H r ) 



9(k T (v,r) - p T )d$ r 



(3.5) 



and k T (v,r) is the transverse momentum of the radiated parton, 



k T (v,r) 



M 2 - 

QQ 

Ax 



(1-X)2 (l-J/2). 



(3.6) 



Equation (3.1) is the analogue of eq. (5.10) of ref. [1]. The function A q (p T ) corresponds to 
A^ LO \p T ) in the notation of ref. [1]. 

4. Generation of the hardest event 

The generation of the hardest event according to eq. (3.1) can be performed in full analogy 
to the case of ref. [6]. We refer the reader to that paper for details. Here we only point out 
the relevant differences with respect to that case. 

The generation of the Born configuration involves in this case two kinds of Born pro- 
cesses, the qq and gg processes. The total cross section is thus given by 



After the Born configuration has been generated, one chooses the process (i.e. g or & given 
flavour of q) with a probability proportional to B g (v,fj, v ), B q (v,fi v ). According to whether 
a q or a g was selected, one follows the same procedure of ref. [6] using to the first or second 
line of eq. (3.1) respectively. 

5. Accuracy of the Sudakov form factor 

Unlike the case of Z pair production, in the case at hand the procedure illustrated in 
section 4 of ref. [6] (i.e. the redefinition of a s given in eq. (4.9) of [6]) is not sufficient 
to guarantee full next-to-leading logarithmic accuracy of the Sudakov form factor. This is 
related to the fact that the heavy flavour production process at the Born level involves more 
than 3 coloured partons [12], so that soft emission cross sections do not simply factorize in 
terms of the Born cross section. Thus, the Sudakov form factor is strictly only accurate to 
leading log. In fact, next-to-leading logarithmic accuracy can be easily recovered at least 
in the large iV c limit, where N c is the number of colours, with a procedure discussed in 
ref. [13]. The implementation of this procedure and the assessment of its impact is left for 
future work. 

6. Colour assignment 

In order to interface POWHEG with a shower Monte Carlo, colour connections must be 
specified. In the case of Z pair production, only one colour structure is present. The 




i 



(4.1) 
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situation is more complex in the case at hand, since more colour structures are relevant. 
This problem is dealt with in exactly the same way as was done in ref. [5], section 6.1. 
We used in fact the same decomposition of the large N c heavy flavour production cross 
section into contributions with different colour structure, and pick the colour structure 
with probability proportional to the value of the corresponding contribution. 

7. Results 

In this section, we present results obtained with the POWHEG method for a choice of 
observables relevant to heavy quark production. Our results will be compared to those 
obtained with MC@NLO. We thus interface POWHEG with HERWIG, in order to make a 
consistent comparison. The formalisms of POWHEG and MC@NLO differ in the treatment 
of contributions of orders higher than NLO, which are beyond the level of accuracy of the 
theoretical computations presently available. The difference is mainly due to the way the 
radiation of matrix element origin is generated, which is typically the hardest radiation. 
Furthermore, scale choices in the two codes are not the same. Therefore, sizable differences 
between the two methods are to be expected for bottom and charm production, where the 
relevant scale is relatively small, while in the case of top production the discrepancies 
should in principle be much less important. However, a detailed comparison between the 
two methods is beyond our present purposes, and it is left for future work. 

We will consider two experimental configurations: pp collisions at y/s = 1960 GeV, 
corresponding to the Tevatron Run II configuration, and pp collisions at y/s = 14 TeV, 
corresponding to the LHC. The results presented in this section have been obtained by 
setting the top and bottom masses equal to 172 GeV and 4.75 GeV respectively. We have 
used the MRST2002 [14] set of parton distributions. When considering the decay of top 
quarks, we have set T t = 1.31 GeV. 

We begin by considering top production. Both the POWHEG and MC@NLO codes 
include the possibility of generating distributions for either undecayed top quarks (which 
we will refer to as "stable top" in the following), or for the decay products of top quarks, 
taking spin correlations into account. We will show examples of both cases. We present in 
fig. 2 the single-inclusive transverse momentum distribution for a stable top quark produced 
at the Tevatron. The ratio between the POWHEG and the MC@NLO results is also shown 
in the lower pane. The POWHEG (solid histogram) and MC@NLO (dashed histogram) 
results are very close to each other over the whole range considered, except in the very 
small p T region, where the POWHEG cross section tends to be larger. 

In fig. 3 we present the distributions of the invariant mass and of the transverse mo- 
mentum of a stable ti pair. The agreement in this case is also quite good; the POWHEG 
p T (tt) distribution is slightly softer than that of MC@NLO. 

The same observables are in an even better agreement in the case of the LHC. This 
is shown in figs. 4 and 5. This can be understood as a consequence of the fact that the 
kinematics of the production process is less constrained than at the Tevatron, therefore 
making the contribution of potentially large logarithms in the perturbative coefficients less 
important. 



- 10 - 





Figure 3: Invariant mass and transverse momentum distributions of tt pairs at the Tevatron. 



We now turn to distributions of the decay products of unstable top quarks. We con- 
sider here only the dilepton channel, which considerably simplifies the analysis. In fig. 6 we 
show two representative single-inclusive distributions, namely the transverse momentum 
and the rapidity of negatively-charged lepton resulting from the t decay, in the Tevatron 
configuration. Again, both the POWHEG and the MC@NLO predictions are shown, to- 
gether with the ratio between the two results. We see that no significant difference is 
present. Similar conclusions hold for the invariant mass and for the transverse momentum 
distributions of the charged-lepton pair, fig. 7, and for the azimuthal distance between the 
two charged leptons, fig. 8. A similarly good agreement is obtained for the same set of 
observables computed in the LHC configuration, figs. 9, 10, and 11. The only visible 
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Figure 4: Transverse momentum distribution of a top quark at the LHC. 




Figure 5: Invariant mass and transverse momentum distributions of tt pairs at the LHC. 



difference is in the overall normalization, which is manifest in figs. 8 and 11. This is due 
to the different choice of scales in the two computations. 

We now turn to the case of bottom production. As is well known, perturbative NLO 
corrections to bottom production are very large, which implies that yet higher-order con- 
tributions are due to play a non-negligible role. As mentioned above, we therefore expect 
that POWHEG and MC@NLO will show larger discrepancies than in the case of top pro- 
duction purely on the basis of fixed-order expansion. There are, however, other sources 
of differences between the two approaches. Although both codes have been interfaced to 
HERWIG in order to obtain the results shown here, the logarithmically-enhanced terms 
beyond the leading one are not the same in the two approaches. Furthermore, if POWHEG 
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Figure 6: Transverse momentum and rapidity of a charged lepton from top decay at the Tevatron. 




Figure 7: Invariant mass and transverse momentum distributions of £ + £ pairs from top decay at 
the Tevatron. 



is interfaced to an MC based on angular-ordered evolution (such as HERWIG), standard 
showers need be supplemented by truncated showers, whose effect is that of restoring colour 
coherence, which is lost because of the requirement that the hardest radiation be always 
the first. Since truncated showers are inherently soft, there are reasons to believe that their 
effects are not too large. At present, the only study of the impact of truncated showers 
has been performed in ref. [15]. There, a POWHEG implementation of e + e~ annihilations 
into hadrons, interfaced to the HERWIG++ Monte Carlo [16], was presented. The effect 
of the truncated shower was found to be small. No studies have been performed in the case 
of hadron collisions. 

In fig. 12 we present sample comparisons between POWHEG and MC@NLO results 
for bottom production at the Tevatron. All observables shown are relevant to lowest-lying 
6-flavoured meson states. We show the single-inclusive p T (upper left pane), the pair p T 
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Figure 8: Distribution of the azimuthal distance between charged leptons from top decay at the 
Tevatron. 




Figure 9: Transverse momentum and rapidity of a charged lepton from top decay at the LHC. 



(upper right pane), and the azimuthal distance, without (lower left pane) and with (lower 
right pane) kinematic cuts; in the latter case, the cuts \y\ < 1 and p T > 5 GeV are 
applied to both the B's of the pair. The two p T distributions show a fair agreement, with 
POWHEG marginally (for single inclusive p T ) or markedly (for the tail of the p T of the 
pair) harder than MC@NLO. There are very clear differences in shape between the two 
azimuthal distributions. The discrepancy tends to be smaller when cuts are applied. As 
for the p T of the pair, POWHEG gives harder results than MC@NLO, which we attribute 
mainly to the different treatment of hard radiation in the two formalisms. Finally, we point 
out that the POWHEG code is capable of producing bottom and charm distributions at the 
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Figure 10: Invariant mass and transverse momentum distributions of £ + £ pairs from top decay 
at the LHC. 
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Figure 11: Distribution of the azimuthal distance between charged leptons from top decay at the 
LHC. 



LHC energy, essentially without negatively- weighted events. In ref. [11] example input files 
are provided for bottom and charm production at the Tevatron and at the LHC, interfaced 
to both HERWIG and PYTHIA. 



8. Conclusions 

In this paper we have presented the implementation of heavy quark pair production accord- 
ing to the POWHEG formalism, which allows an NLO QCD computation to be matched 
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Figure 12: Bottom hadron distributions at the Tevatron. 

with Parton Shower simulations. The fortran code we have constructed can be used to 
predict any infrared safe observable in ti, bb, and cc production at hadron colliders. We 
have compared our results with MC@NLO for ti production at the LHC and at the Teva- 
tron, and for bb production at the Tevatron. In the case of top production, we observe a 
very good agreement between POWHEG and MC@NLO, for all the observables we have 
considered. On the other hand, the two approaches differ significantly for some observables 
in bb production, which implies that for such low-mass quarks perturbative corrections of 
order higher than next to leading are likely to play a non-negligible role. 

In general, the agreement of the MC@NLO and POWHEG approaches is quite re- 
markable, in view of the fact that the two methods differ considerably in several aspects, 
summarized below: 

1. The Sudakov form factors in MC@NLO and POWHEG are different: MC@NLO uses 
HERWIG's Sudakov form factor, POWHEG has its own (see eqs. (3.4) and (3.5)). 

2. The hardest emission in POWHEG carries a strong coupling evaluated at the p T of 
the emission, and a Sudakov form factor. On the other hand, in MC@NLO only S 
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events, 2 that have all emissions entirely performed by the shower, have these features. 
The H events are evaluated at the scale of the hard process, and they carry no Sudakov 
damping for small transverse momenta. This difference may show up for relatively 
low transverse momentum. We remind the reader that EI events can have negative 
weight, so that it is difficult to understand in which direction this difference affects 
the results. 

3. The POWHEG approach lacks the truncated showers. 

4. Subleading terms in the shower may differ in the two approaches, due to the reshuf- 
fling of the splitting processes in the shower illustrated in ref . [1] . 

5. The EI events in MC@NLO may be followed by radiation, generated by the shower, 
with a p T harder than the p T of the EI event. In POWHEG harder emissions from 
the shower are always vetoed. 

Because of the many differences, it is also difficult at this stage to understand what causes 
the differences in the distributions we have presented. Here we just make a few speculations 
about the possible origin of the differences, that should only be taken as hints for further 
studies. First, we look at top production. We see there that the inclusive p T spectrum, 
the mass of the pair and the p T of the pair differ in the very small p T or m(tt) region, 
POWHEG being generally higher (see figs. 2 to 5). The lack of soft-truncated showers 
in POWHEG (item 3 of the above list) could possibly cause this effect. On the other 
hand, soft radiation is also treated differently as far as the hardest emission is concerned, 
as specified in item 2 of the above list. The fact that the difference goes in the opposite 
way (i.e. POWHEG is below MC@NLO) for the transverse momentum of the top pair 
at LHC also shows that the lack of truncated showers (that would lower the POWHEG 
distribution) cannot be the whole answer. A second effect we notice is the considerable 
difference in the azimuthal distance of the bottom pair (see fig. 12) especially in the region 
where the two heavy mesons are near in azimuth. Here, the lack of soft-truncated showers 
in POWHEG, or differences in subleading shower effects in the two methods, could yield a 
different degree of smearing of the azimuthal distance. The POWHEG result seems to have 
more smearing than the MC@NLO result. The lack of soft-truncated showers is likely to 
have the opposite effect. Also, the faster rise for small azimuthal difference is unlikely to be 
due to shower effects in general, and would suggest to look for effects in the hard radiation 
mechanism. Thus, items 2, 4 and 5 may be responsible for these differences. The third 
effect we consider is the harder tail of the p T distribution of bottom pairs in POWHEG. 
This can only be ascribed to genuine higher order effects in the hard emission, as may arise 
from items 1 and 2. 3 

2 See refs. [4, 5] for the definition of § and H events. 

3 In a yet unpublished revision of the MCQNLO code the p T spectra for B production turn out to be 
harder than the ones shown here. The azimuthal distance distributions do not differ from those presented 
here. 
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